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The Bose-Einstein (or HBT) correlation functions are evaluated for the 
fractal structure of QCD jets. These correlation functions have a stretched 
exponential (or Levy-stable) form. The anomalous dimension of QCD de- 
termines the Levy index of stability, thus the running coupling constant 
of QCD becomes measurable with the help of two-particle Bose-Einstein 
correlation functions. These considerations are tested on NA22 and UAl 
two-pion correlation data. 

PACS numbers: 05.40.Fb, 13.85.Hd, 13.87.Ce, 25.75.-q, 25.75.Gz 



1. Introduction 

The study of fractal phenomena was introduced to high energy particle 
and nuclear physics by Bialas and Peschanski in ref. PP, see also ref.|2] 
for a review. In QCD, jets emit jets that emit additional jets and so on. 
The resulting fractal structure of QCD jets was explored with the help of 
a beautiful geometric picture in refs. ^|1J|5]. These ideas were developed 
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further by the Lund group in refs. [HIEIIHI; by Dokshitzer and Dremin in 
ref. IHl as weU as by Ochs and Woisek in refs. |lUllllj . Both theoretical and 
experimental aspects of the so-called intermittency or fractal structures in 
high energy physics were reviewed by De Wolf, Dremin and Kittel in ref. jl2j . 

Bialas realized, that Bose-Einstein correlations and intermittency might 
be deeply connected US]- The mathematical properties of Bose-Einstein 
correlation functions for Levy stable (convolution invariant) sources were 
written up by three of us in refs. jl41 115j . Here we add a physical interpre- 
tation and show, that the fractal properties of QCD cascades can naturally 
be measured by the Levy index of stability of Bose-Einstein correlations. 
Our analytical results are similar in spirit to the numerical investigations of 
Wilk and collaborators in ref. ^16) . 



2. (Multi) fractal structure of the QCD jets 

In this section we recapitulate earlier theoretical results of the Lund 
group, 131 13 El El 13 El ) that related the properties of QCD cascades to 
intermittency. These results are based on a beautiful geometric interpreta- 
tion of the color dipole picture, and on an infrared stable measure on parton 
states related to hadronic multinlicitv. 




Fig. 1. The phase-space of QCD jets in the {y, k) plane, where k = log(A:( ). (a) The 
phase-space available for a gluon emitted by a high energy qq system is a triangular 
region in the {y,K) plane, (b) If one gluon is emitted at the phase-space 

for a second (softer) gluon is given by the area of this folded surface, (c) The total 
gluonic phase-space can be described by this multifaceted fractal surface O^E]. 

A high energy qq system radiates gluons according to the dipole formula 
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hence the phase-space for the emission of a gluon is given by the relation 

(2) 
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which corresponds to the triangular region in a {y,\nk\) diagram as shown 
in Fig. 1 (a). If two gluons are emitted, then the distribution of the hardest 
gluon is described by eq. (1). The distribution of the second, softer, gluon 
corresponds to two dipoles, the first is stretched between the quark and the 
first gluon, and the second between the first gluon and the anti-quark. The 
phase-space available for the second gluon corresponds to the folded surface 
in Fig. 1 (b), with the constraint A;^ 2 

j^, as the first gluon is assumed to 
be the hardest one. This procedure can be generalized so that the emission 
of a third, still softer gluon corresponds to radiation from three color dipoles, 
with n gluons emitted already the emission of the (n -|- l)-th gluon is given 
by a chain of n -|- 1 dipoles. Thus, with many gluons, the gluonic phase- 
space can be represented by a multi-faceted surface as illustrated in Fig. 1 
(c). Each gluon adds a fold to the surface, which increases the phase-space 
for softer gluons. (Note, that in this process the recoils are neglected, as 
is normal in leading log approximation). Due to its iterative nature, the 
process generates a Koch-type fractal curve at the base-line. The length 
of this base-line of the partonic structure on Figure l.c is proportional to 
the particle multiplicity. This curve is longer, when studied with higher 
resolution: it is a fractal curve, embedded into the four-dimensional energy- 
momentum space, characterized by the fractal dimension 

= l + (3) 

or one plus the anomalous dimension of QCD |3J|2|S]. 

With the help of the Lund string fragmentation picture, this fractal in 
momentum space is mapped into a fractal in coordinate space, and the 
constant of conversion is the hadronic string tension, k ~ 1 GeV/fm. This 
mapping does not change the fractal properties of the curve. The emission 
of softer and softer gluons corresponds to a smaller and smaller modification 
of this curve, as a gluon with a very small transverse mass creates a very 
small kink on the Lund string. Hence this process is infrared stable. 

3. Bose-Einstein correlations for Levy stable source distributions 

Let us discuss here the stability of the particle emitting source in QCD, 
and consider the Bose-Einstein correlation functions for such sources. The 
two-particle Bose-Einstein correlation function is defined as the ratio of the 
two-particle invariant momentum distribution to the product of the single- 
particle invariant momentum distributions: 



(4) 
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If long-range correlations can be neglected or corrected for, and if the 
short-range correlations are dominated by Bose-Einstein correlations, this 
two-particle Bose-Einstein correlation function is related to the Fourier- 
transformed source distribution. For clarity, let us consider the case of a 
one-dimensional, factorized coordinate and momentum space distribution, 

S{x,k)=fix)g{k). (5) 

In this case jl41 115j , the Bose-Einstein correlation function is 

C2{h,k2) = l + \f{q)\'', (6) 

where the Fourier transformed source density (often referred to as the char- 
acteristic function) and the relative momentum are defined as 

/(,) = /d.exp(i,.,/(.), , = fe-te. (7) 

Let us focus on the property of particle emission from QCD jets, that 
the fractal defining the particle emission is infrared stable: adding one more, 
very soft gluon does not change the resulting source distributions. Thus the 
source of particles is stable for convolution. The Bose-Einstein correlation 
functions for such particle emitting sources were evaluated recently by three 
of us, which we summarize in this section following refs. Jl] and |15j . 

For the case of the jets decaying to jets to jets and so on, the final 
position of a particle emission is given by a large number of position shifts, 
hence the distribution of the final position x is obtained as a convolution, 




Various forms of the Central Limit Theorem state, that under certain con- 
ditions, the distribution of the sum of a large number of random variables 
converges (for n oo) to a limit distribution. In case of "normal" elemen- 
tary processes, that have finite means and variances, the limit distribution 
of their sum is a Gaussian. Stable distributions are precisely those limit 
distributions that can occur in Generalized Central Limit theorems. Their 
study was begun by the French mathematician Paul Levy in the 1920's. 
The stable distributions are frequently given in terms of their characteris- 
tic functions, as the Fourier transform of a convolution is a product of the 
Fourier-transforms, 

n 
i=l 
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and limit distributions appear when the convolution of one more elementary 
process does not change the form of the limit distribution, but it results only 
in a modification of its location and scale parameters. The characteristic 
function of univariate and symmetric stable distributions is 

f{q)=e^v{iq5-\iqn, (10) 

where the support of the density function /(x) is (—00,00). Deep mathe- 
matical results imply that the index of stability, a, satisfies the inequality 
< a < 2, so that the source distribution be always positive. These Levy 
distributions are indeed stable for convolutions, in the following sense: 



fi{q) = exp {iq6i - iTigD , J| fi{q) = exp {iq5 - |7g|'*) , (11) 

n n 

7" = E^f, S = Y.5,. (12) 

1=1 i=l 

Thus the Bose-Einstein correlation functions for univariate, symmetric 
stable distributions (after a core-halo correction, and a re-scaling) read as 

C(g;a) = 1 + Aexp(-|gi?|"). (13) 

For the special value of a = 2 we obtain the well known Gaussian case. 
Refs. ^1] and ^3] discuss further examples and details and generalize these 
results to three dimensional, hydrodynamically expanding, core-halo type 
sources as 

C{q,K;a) = 1 + A(K) exp |^-| ^ 9^4 W^il"/' j • (14) 

Figure 111 (a) shows univariate symmetric Levy source distributions. Fig- 
ure [21 (b) illustrates, that values of the index of stability a are related to 
the tails of these distributions. If a < 2, for large values of s = r/R the 
Levy sources decay as fa{s) oc Figure|21(c) shows a two-dimensional 

symmetric Levy stable source. Bose-Einstein correlation functions for Levy 
stable source distributions are shown in Figure El (d) for various values of 
the index of stability a, and for a constant value of the radius parameter 
R. These Bose-Einstein correlation functions are sensitive to the value of a 
not only in the small Q < fi/R region, but are also in the "large" relative 
momentum region of Q > %/ R. Thus these correlations are sensitive to the 
structure of the particle emission in the region which is shaped by the jets. 

A random walk, where the length of the steps is given by a Levy distri- 
bution, and the direction of the steps is random, corresponds to a fractal 
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Log[f(s)] Levy sources, lin-log scale 




Fig. 2. (a - top left) Source functions for univariate symmetric Levy laws, as a 
fmiction of the dimensionless variable s = r/R, on a linear- linear scale, for various 
values of the Levy index of stability, a. (b - top right) Same as (a), but on a linear- 
logarithmic plot. Note the power-law tails that decay in case of a < 2 for large 
values of s as /a(s) oc s^^^". (c - bottom left) Source function for two-dimensional 
symmetric Levy laws, for a — 0.6. (d - bottom right) Bose-Einstein correlation (or 
HBT) correlation functions for univariate symmetric Levy laws, for a fixed scale 
parameter of ii = 0.8 fm and various values of the Levy index of stability, a. 



curve, in physical terms it can be interpreted as the path of a test particle 
performing a generalized Brownian motion. This motion is referred to as 
anomalous diffusion and the probability that the test particle diffuses to 
distances r greater than a certain value of |s| is given by P{r > \s\) oc |s|~"- 
This relation is valid for anomalous diffusion not only in one, but also in 
two and three dimensions. Thus the Levy index of stability a is the fractal 
dimension of the trajectory of the corresponding anomalous diffusion jl7j . 
When we apply this result to QCD, there are two key considerations. 

First, if gluon radiation is neglected, the qq system hadronizes as a 
dimensional hadronic string, which has no fractal structure. If the gluon 
emission is switched on, the emission of gluon n from one of the n dipoles 
corresponds to a step of an anomalous diffusion in the plane transverse to 
the given dipole. Hence the anomalous dimension of QCD equals to the 
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Levy index of stability of this anomalous diffusion, 

y^="Levy (15) 

Second, data on Bose-Einstein correlations are often determined in terms 
of the invariant momentum difference Qmv = \/— (pi — P2)''^- Bose-Einstein 
correlation functions that depend on this invariant momentum difference 
can be obtained within the framework of the so-called r-model. This model 
assumes a broad proper-time distribution, H{t) and very strong correlations 
between coordinate and momentum in all directions, x^/t oc p^/m(^t) ■ Hence 
(xi —X2)ipi —P2) oc rQfj^^, see refs. ^lEl for details. In this case, the Bose- 
Einstein correlation function measures the Fourier-transformed proper-time 
distribution H in the following, unusual manner: 

C2iQin.)^l + XReH^^], (16) 

where m(() stands for the (transverse) mass of the pair for (two)- or more 
jet events. From this relation it follows, that Obec = '^ctLcvy Thus we 
find the following relationship between the strong coupling constant and 
the exponent of an invariant relative momentum dependent Bose-Einstein 
correlation function: ^ 

as = goLc- (17) 

4. Application to NA22 and UAl data 

In ref. [Tij three of us have fitted the NA22 pD] and the UAl data [H] on 
two-particle Bose-Einstein correlation functions with the Levy stable form 
of eq. (|13l) . The results were summarized in Table 1 of that paper. Here 
we re-interpret the exponent of this fit with the help of eq. (|17|) and extract 
Us, the coupling constant of QCD, as given by Table 1 of this paper. 



Parameter 


NA22 


UAl 


Value Error 


Value Error 


Obec from ref. [lij 


0.67 ± 0.07 


0.49 ± 0.01 


as from eq. (flTl) 


0.24 ± 0.05 


0.13 ± 0.01 



Table 1. Converting the best Levy fits to UAl and NA22 two-particle correlations 
using the r-model converted with ea. l|17() to values for the running QCD coupling 
constant as- 
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5. Summary, conclusions 

Using the picture of strongly correlated coordinate and momentum space 
distributions, we determined the (running of the) strong coupling constant 
from NA22 and UAl two-pion Bose-Einstein correlation measurements. 
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